Conservation and persistence of spin currents and their relation to the 
Lieb-Schulz-Mattis twist operators 
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Systems with spin-orbit coupling do not conserve "bare" spin current j. A recent proposal for a 
conserved spin current J [J. Shi et.al Phys. Rev. Lett. 96, 076604 (2006)] does not flow persistently 
in equilibrium. We suggest another conserved spin current J that may flow persistently in equilib- 
rium. We give two arguments for the instability of persistent current of the form J: one based on 
the equations of motions and another based on a variational construction using Lieb-Schulz-Mattis 
twist operators. In the absence of spin-orbit coupling, the three forms of spin current coincide. 

PACS numbers: 72.10.Bg,72.20.Dp, 73.63.Hs 



I. INTRODUCTION 

Spintronics aims to manipulate the electron's spin in- 
stead of its chargeJ^i^ii A requisite for spintronic de- 
vices is the generation of a spin current. ^'^ In recent 
years, exciting new proposals and predictions have been 
made regarding the possibility of dissipationlcss spin 
currents- primarily through the use of the Spin Hall Ef- 
fect (SHE) J^iii The SHE may generate sizable HaU spin 
currents and open the door for the realization of many 
of these notions. Observation of the SHE have been re- 
ported by employing the Kerr effect in n-doped GaAs-- 
and a p-n junction light emitting diode systemfi^, Much 
fundamental work on spin currents has been carried out, 
e.g., Ref . iiiHilil. 

A fundamental problem in spintronics is the non- 
conservation of the bare spin currents (the sum of the 
product of each electron spin multiplied by its velocity, 
see Eq. (|31))- This non-conservation is triggered by spin- 
orbit effects. Charge conservation and the continuity 
equations that it gives rise to are of crucial importance 
in electronic circuits. Is there an analogue for spins? The 
quest for conserved spin currents has been the subject of 
much attentioUfiSiii^i^ An early step was made in Ref.fll 
that proposed a conserved spin current by augmenting 
the bare spin currents by an additional contribution so 
as to include the effects of torques triggered by spin-orbit 
effects. The work of Ref. [H assumed that the average 
torque vanishes in the bulk and that the "torque dipole 
density" on the system's surface can also be neglected. 
The "torque dipole density" was a field whose divergence 
led to the actual torque acting on spins at a particular 
location. This proposal was exciting. The spin accumu- 
lation observations of Ref. [Ill and Ref. [H indeed cannot 
measure the bare spin current. 

We suggest that the notions of conservation and of per- 
sistence are independent for spin currents. For example, 
the bare spin current in Eq. ([3]) in general is not con- 
served, and may persist in equilibrium (See Eq. (|10p ). 
Rashba noticed the possibility of persistent equilibrium 
flow of the bare spin current for a particular class of 



spin-orbit coupled modelsJi Here we present a model- 
independent argument for the result. In contrast to the 
bare spin current, the conserved spin current J given by 
Eq. (12) does not flow in energy eigenstates so any persis- 
tent flow implies the system is not in equilibrium. Fur- 
thermore, we have found another conserved spin current 
J and associated spin density p that may flow persis- 
tently. When spin is a good quantum number, these dis- 
tinctions collapse: the bare spin current is conserved, and 
none of the three spin currents flow persistently. How- 
ever, in a generic, spin-orbit coupled system, the defini- 
tions of spin current differ with respect to their conser- 
vation and persistence. 

In Section (jlll, we recall the standard definition of spin 
density and bare spin current. We then review the con- 
served spin current proposed in Ref. [3, and contrast it 
with an alternative, conserved spin current (See Eq. (O). 
Section (IIIip further distinguishes this novel, conserved 
current by computing its expectation in a general spin- 
orbit coupled system. We show that, in contrast to the 
previous proposal, it may fiow persistently, just like the 
bare spin current. In Section pvp we give a physical ar- 
gument to explain this distinction: any state that carries 
the conserved spin current proposed by Ref. 1^ may dissi- 
pate energy by relaxing to a lower energy state, which we 
explicitly construct. The argument extends the one used 
by Bohm to establish that persistent charge current flows 
only in meta-stable states, not the ground-statci^ Fur- 
ther, the construction generalizes one used by Lieb et al. 
to study quasi-one-dimensional, spin-1/2 systemsi^^ We 
conclude with a summary and discussion of our results. 



II. DEFINITIONS OF SPIN CURRENT AND 
SPIN DENSITY 

Let us consider a general, interacting system with 
Hamiltonian i7[{ri}, {pi}, {ctq}]. Here, {ri} and {pi} de- 
note the phase space variables and {ctq} reflect the spin 
states of all non-scalar particles in the system. The spin 
density along a given quantization axis (taken to be the 



2 



-axis) has the form 



(1) 



where, is the particle density, and erf is the spin oper- 
ator along the 2;— axis, given by a Pauli matrix for spin 
one-half particles. In Reffisl Shi et. al. suggested the 
following form for the spin current: 



(2) 
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III. CONSTRAINTS ON SPIN CURRENT 
FROM EQUATIONS OF MOTION 

As noted in Rcf. Il8|, the form of the spin current in 
Eq. 12]) does not flow in energy eigenstates. We give a 
simplified derivation, then extend it to finite tempera- 
ture, and use the result to deduce constraints on equilib- 
rium flow of the other two forms of spin current Eq. ([3|) 
and Eq. 

The spin displacement operator 



L = ^ XiCrf 



(7) 



where, the spin torque is given by = dcrf /dt — 
1/ {ih)[af ,H]. In the absence of spin-orbit coupling, the 
spin torque vanishes, and the spin current J coincides 
with the "bare" spin current 



(3) 



dxi /dt 



where, the velocity is given by 

In contrast to the bare spin current the current 
in ([2|) satisfies the continuity equation, even in systems 
with spin-orbit couplingjiS 



(4) 



Physically this means that in a system with boundary, 
the flow of spin current ([2]) across the boundary into the 
system leads to spin accumulation within the system. We 
note that conservation of J only occurs when the average 
torque vanishes: J2i{'''i) — ^-^^ 

Interestingly, it is possible to construct another spin 
current J and spin density p that satisfy the continuity 
equation. Consider the following spin current: 
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(5) 



where, the omitted terms are higher order in the spin- 
orbit coupling. We define a corresponding spin density 
p = erf (i = 0)ni by evaluating at the same reference 
time. Remarkably, this form of the spin density and spin 
current also satisfies the continuity equation: 



9p 

dt 



V-J = 0. 



(6) 



To see this, notice that the label crf(t = 0) does not 
evolve, but merely labels each particle. From the con- 
tinuity of charge, the continuity of these labels follows. 
We remark that the current J obeys the continuity equa- 
tion ([6]) even when the average torque does not vanish, in 
contrast to the current J. In the appendix, we perform 
the expansion in spin-orbit coupling explicitly. 



allows us to evaluate the spin current in Eq. ([2]): 
{n\3\n) = {n\^\n) = -t{n\[L,H]\n) = 0. 



(8) 



Here, \n) denotes the nth eigenstate of H with energy En- 
The second equality in Eq. ([8]) follows from the equations 
of motion while the third follows as (n| and |n) refer to 
the same eigenstate of H with specific eigenvalue (£"„). 
Here and henceforth, we set h = 1. 

As an immediate corollary, we notice that the spin cur- 
rent in Eq. ([2]) vanishes at finite temperatures as well: 



(J) 



2„(n|J|n)e- 



-fJE„ 



0, 



(9) 



where, the inverse temperature is (3 — l/ksT. 

If wc now invoke Eqs.lO [9]), we have the following con- 
straint on the bare spin current Eq. ([3|) in equilibrium: 



(10) 



This is so due to Eq. ([9]) , and the relation bewteen Eq. ([2|) 
and Eq.(I3]). Similarly, the conserved current J may have 
a non-vanishing expectation value even at equilibrium: 



(J) = -2E(r^x,) + 



(11) 



where, the omitted terms are higher order in the spin- 
orbit coupling. When spin is a good quantum number, 
all three forms of the spin current vanish in equilibrium. 
Only in the presence of spin-orbit coupling do the equi- 
librium currents differ. 



IV. A VARIATIONAL PROOF USING THE 
LIEB-SCHULZ-MATTIS TWIST OPERATOR 

We now regress and derive a version of Eq. ([8]) for the 
ground state sector by a generalization of the variational 
arguments of Ref. [2l|. This proof makes contact with 
topological Lieb-Schulz-Mattis^^ type twist operators. 
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We consider the ground state |0) of a general system 
with spin-orbit coupHng and prove by a variational argu- 
ment that within it the spin current vanishes: 

(0|J|0) = 0, (12) 

where, we take spin current J as in Eq.Q. By definition, 
the ground-state energy 

Eo = (0|i/|0). (13) 

is the lowest attainable energy and lies at the bottom of 
the spectrum. Let us assume that we have a state lip) 
with a finite (i.e. non-zero) spin current. We will then 
prove that cannot be the ground-state, by construct- 
ing another state \ip{6P)) with a lower energy. That is, 
if (■0|J|'0) ^ then, the energy E{SP) of the "twisted" 
state \ip{S'P)) can be made to be lower than the energy 
E{0) of by choosing the variational parameter SP 
appropriately: 

E{5P) < E{0). (14) 

In this way we establish that the ground-state |0} cannot 
carry a spin current of the form in Eq. ([2]), since, by 
definition, it has the lowest energy. 

We construct the variational state \ip{SP)) by applying 
a spin-dependent boost to the state 1-0): 

\^) = exp(-i5P-L)|?A) 

= exp(-z(5P.^x,an, (15) 

i 

where, we have used the spin-displacement operator L 
from Eq. ([7]). Compared to the momentum of each 
particle in |?/'((5P)) with spin up = +S changes by 
-|-S'(5P while the momentum of particles with spin down 
erf — —1 changes by — ^^P. The momentum boost SP 
will be our variational parameter. Lieb-Shulz-Mattis per- 
form a similar construction but restrict the system geom- 
etry to be quasi-one-dimensional, and the spin S to be 
half-integer.^^ We require no such restrictions: the con- 
clusion Eq. (fT2)l holds in one-, two- and three-dimensions, 
for systems with half-integer spin particles as well as 
those containing integer spin particles. 

The energy of \'ip{SP)) for small dP may be expanded 

{^{SP)\H\i:{SP)) = imii^) ~ tSP ■ {nH,J2a',^iM 

i 

+ 0{{5Pf). (16) 

Now invoking the equation of motion (dL/dt) — i{[H, L]) 
and the definition of the spin current J — dL/dt in Eq. ([2|) 
we have 

E{dP) = E{0) + 5P ■ (V^|J|V^) + 0{{5Pf). (17) 

It follows that if the state \^) carries a spin current 
(V-'|J|'0) 7^ then there always exists a sufficiently small 
5P such that E{5P) can be made smaller than i?(0), by 
choosing the direction of the momentum boost for spin- 
up particles 5P opposite to the current direction (J). 
Since the ground-state has the lowest energy, we must 
have Eq. a special case of the more general Eq.®. 



V. CONCLUSIONS AND DISCUSSION 

In spin-orbit coupled systems, the definition of spin 
current has inherent ambiguityi^ By redefining the cur- 
rent, it is possible to change its properties. We consider 
three definitions of the spin current that generically differ 
with respect to the two properties that interest us most, 
continuity and persistence. The bare spin current j de- 
fined in Eq. ^ may persist in equilibrium, but generally 
does not lead to spin accumulation. The spin current J 
defined in ref. il8j obeys the conservation law ([4]) when 
the average torque vanishes. However, the expectation 
{n\3\n) vanishes in all energy eigenstates \n) of a system 
with a time-independent Hamiltonian. A third definition 
([5]) of the spin current J, obeys a continuity equation with 
a modified spin density p, and may flow persistently in 
equilibrium. 

Expanding in powers of the spin-orbit coupling, we find 
that J differs at leading order from both the bare spin 
current j and the conserved current J of Ref. [ll|. This 
means that in principle, the expectation value (n|J|n) 
need not vanish, and, a conserved spin current J may 
flow persistently even in equilibrium despite spin-orbit 
coupling. We conclude that in general spin-orbit coupled 
systems, no definite relation exists among the concepts 
of persistence and conservation of spin current. 

We give a physical argument for the instability of per- 
sistent spin current of the form J defined in ref. 18. The 
argument makes contact with the Bloch's famous result 
for the instability of persistent charge currents. In pass- 
ing, we note a fundamental relatioiij_summarized in (jlSp . 
between J the spin current of Ref. [l8' and the twist gen- 
erator used by Lieb, Schulz, and Mattis. It would be 
interesting to formulate this relation in an SU (2) invari- 
ant fashion. 

Several experiments bear directly on these 
issueSi ^^i^^'^^ Transport measurements in HgTe quantum 
wells provide evidence for gapless edge modes that are 
protected from non-magnetic impurity scattering and 
that respond singularly to an external magnetic fieldi^^ 
This is consistent with a "helical" picture for the edge 
modes, with a bare spin current j accompanying the 
charge current. Spectroscopic measurements of the 
surface excitations in Bii-^^Sb^; find an odd number of 
gapless Dirac points on the surface, suggesting topo- 
logical stability of the gapless modes to non-magnetic 
impurities.— Further, in GaAs/GaAlAs quantum wells, 
time-resolved Kerr rotation experiments see evidence for 
a long-lived excitation mode that transports both spin 
and charge.^- 

These systems contain strong spin-orbit coupling, so 
the bare spin current j accompanying the gapless charge 
modes need not give rise to a spin accumulation. A de- 
tectable change in spin density p, could arise in response 
to an applied bias due to the flow of the current J, pro- 
vided that the average torque density vanishes. For pure 
Rashba spin-orbit coupling, non-magnetic impurity scat- 
tering eliminates spin accumulation even under an ap- 
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plied bias.-^ The sy terns described above are stable to 
non-magnetic impurities, opening up the possibility of 
detectable spin accumulation. 

We gratefully acknowledge discussions with A. V. Bal- 
atsky. ZN thanks the CMI of WU for support. NBA 
is grateful to Los Alamos National Laboratory and the 
National Center for Theoretical Sciences for their hospi- 
tality. 

APPENDIX A: A DERIVATION OF CONSERVED 
SPIN CURRENT J 

We begin by considering the conservation of charge 
current. In the absence of spin-orbit coupling, this case 
also describes spin. Charge conservation and gauge in- 
variance are tightly linked. To see this in the simplest 
context, we consider a system with approximate d + 1- 
dimensional Lorentz invariance such as graphene (d — 2) 
or bismuth (d = 3), and adopt the relativistic notation 
A^, where, Aq = (p is the scalar potential and Aa for 
a = 1, 2, ... d are the components of the vector potential. 
For convenience, we adopt the units h = c = 1. 

Let us now suppose that the low-energy effective action 
SefflAfj,] is invariant under a transformation of the gauge 
field 

Af,^Af, + SAf, ^ Af, - a^A (Al) 

with A(a;p) any smooth function. To see how this con- 
strains charge flow, we formally expand the low-energy 
effective action in powers of the gauge field A^,, 

Sef f [A^] = J d'^xdtJc^'A^ + ..., (A2) 

where, the charge current is given formally by the func- 
tional derivative — ^^sI^^jA, By definition of the 
charge current, the gauge transformation (jAip leads to 
the following change in the effective action: 

SSeff = J d'^xdtJc^dA^, 

= J d'^xdt{df,J^)A+ (A3) 

where, we have integrated by parts in the second line 
and dropped the contribution from the surface of the 
system. Gauge invariance requires SS^ff = for arbi- 
trary, smooth A, which can only be satisfied if the charge 
current is conserved: 

d^,J^ = 0. (A4) 

Conversely, charge conservation (jA4p implies invariance 
of the effective action under the gauge transformation 
(|A1[) . Thus, charge conservation and gauge invariance 
have an intimate relationship. 

To make contact with our results in Section II, let us 
express this connection between gauge invariance and 



charge conservation in the language of operators. For 
a system with time independent charges (e^) the charge 
density operator is given by 

Pc^^eiUi, (A5) 

i 

with rii the number density operator. If we now assume 
gauge invariance and the usual Maxwell dynamics for 
the gauge field, then we find charge is conserved at the 
operator level, 

% + V • Jc = 0, (A6) 

ot 

with the charge current operator given as follows: 

i 

Here is the velocity of the i-th charge carrier in the 
system. 

For example, consider a single-particle with mass m 
and charge e in state lip). Symmetrizing (|A7p . we find 
the following familiar expression for the charge current 
(Here a = 1, 2, ... d run over the spatial components): 

{i;\j^\ij) = e I d''x[r[^^Da]^-^[^^D,]r)m 

with the spatial components of the covariant derivative 
given by, Da = -^^—ieAa- The velocity operator is given 
by V = D/m, and we see here that the operator expres- 
sion (jA7p describes the usual charge current in (jASp . 

Let us now apply this to spin systems. As the spin 
states change with time, we obviously cannot simply sub- 
stitute Ci <^i{i) in order to get a conserved current. If 
we set, however, — > tTf(i = 0) then the steps above 
can be reproduced word for word. Although, we will re- 
fer to spin currents in what follows, our results hold for 
any other "charge" field a which is allowed to vary in 
time. This is so as cr^(O)- the value of the z component 
of the spin at time t — Q- is, by definition, time indepen- 
dent. We may relate cr^(O) to a^'{t) by a time reversed 
evolution, 

a-(0) = e-*/o a^(t)e^/o (A9) 

Here, e is a time ordered exponential. The to- 

tal number operator rii (the sum over both up and down 
spin flavors) and its associated number current satisfy 
the continuity equation (|A7[) . so the continuity equation 
for the spin current ([5]) given by ([5]) follows. 

It is interesting to formally expand the spin density p 
in powers of the spin-orbit coupling: 

i 

i i 
i i 
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To zeroth order, p coincides with the bare spin density 
p. However, already at leading order the two densities 
differ. 

Turning now to the spin current, we can use the ex- 
pression in (|A7[) to write down a corresponding conserved 
spin current. Again it is interesting to formally expand 
the expression in powers of the spin-orbit coupling: 

i 



i i 

+ 7E^v.--^^v. + .... (All) 

i i 

We find that to zeroth order in the spin-orbit coupling, 
J coincides with the bare spin current j in ([3]) and the 
conserved spin current J in ([2]). At leading and higher 
order, however, the three currents differ. 
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